Charge localization and phonon spectra in hole doped La 2 NiC>4 
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The in-plane oxygen vibrations in La2NiC>4 are investigated for several hole-doping concentrations 
both theoretically and experimentally via inelastic neutron scattering. Using an inhomogeneous 
Hartree-Fock plus RPA numerical method in a two-dimensional Peierls-Hubbard model, it is found 
that the doping induces stripe ordering of localized charges, and that the strong electron-lattice 
coupling causes the in-plane oxygen modes to split into two subbands. This result agrees with 
the phonon band splitting observed by inelastic neutron scattering in L^-iSr^NiCU. Predictions 
of strong electron-lattice coupling in La^NiCU, the proximity of both oxygen-centered and nickel- 
centered charge ordering, and the relation between charged stripe ordering and the splitting of the 
in-plane phonon band upon doping are emphasized. 

PACS numbers: 75.60.Ch, 74.25.Kc, 71.45.Lr, 71.38.+i 



There is currently great interest in the importance of 
charge localization and ordering tendencies in a variety 
of doped transition metal oxides: including nickelates, 
bismuthates, cuprates, and manganites [0-0. Recent ex- 
periments have suggested nanoscale coexistence of charge 
and spin ordering, as well as related multiscale dynamics 
p]-^0|. The cuprates have been widely investigated, both 
theoretically and experimentally, as this inhomogeneity 
may be related to high-temperature superconductivity 

The nickelates are considered strong electron-lattice (e- 
1) coupling systems, which helps stabilize charge-ordering 
in the form of " stripe" phases . For the commen- 

surate 1/3 doping case of Lai.67Sro.33Ni04, it has been 
shown in optical absorption and Raman scattering exper- 
iments that new phonon modes appear when the tem- 
perature is lowered below the stripe-ordering tempera- 
ture (T SO —240 K); this is a signature of the stripe forma- 
tion p-[To||. Until now, only the temperature dependence 
and the apical oxygen (Ni-0(2)) vibrations have been 
investigated. However, the doping dependence and the 
in-plane oxygen vibrations (Ni-O(l) stretching modes) 
are also very important for the properties of the quasi- 
two dimensional nickelate materials. In this study, we 
use an inhomogeneous Hartree-Fock (HF) plus random- 
phase approximation (RPA) numerical method for a two- 
dimensional (2D), four-band Peierls-Hubbard model, to 
interpret the inelastic neutron scattering spectra. This 
reveals specific signatures of the stripe patterns in the 
in-plane oxygen phonons. 

Our main results are: (i) There is agreement between 
the results from our multiband model including electron- 
electron and e-1 interactions and the inelastic neutron 
scattering spectra for the in-plane oxygen vibrations with 
various commensurate hole-doping concentrations; (ii) 
The theoretical results predict new vibrational modes 



("edge modes") which are associated with oxygen mo- 
tions near localized holes or in the vicinity of stripes; 
(iii) The e-1 coupling strength, at which the best agree- 
ment between our model and the inelastic neutron scat- 
tering data is achieved, is close to the transition from 
an oxygen-centered stripe phase to a nickel-centered one. 
This suggests that the nickelates may be in a mixed state 
of both stripe phases, and sensitive to temperature, pres- 
sure and magnetic field. 

The inelastic neutron scattering spectra were measured 
on polycrystalline La2- 2; Sr 2 ;Ni04 for various doping con- 
centrations, x = 0, 1/8, 1/4, 1/3, 1/2. Time-of-fhght 
neutron scattering measurements were performed on the 
Low Resolution Medium Energy Chopper Spectrometer 
at Argonne National Laboratory's Intense Pulsed Neu- 
tron Source. For all measurements, an incident neutron 
energy of 120 meV was chosen and data were summed 
over all scattering angles from 2° — 120°. Detailed in- 
formation on the experiment, as well as the preparation 
of the samples used, can be found in Ref. |21|]. The fo- 
cus is on the (neutron-scattering-weighted) generalized 
density-of-states (GDOS) of phonons and particular at- 
tention is given to the in-plane oxygen vibrations, i.e., 
the Ni-0(1) stretching modes. 

Figure [l] shows the experimental GDOS for several hole 
concentrations at T=10 K for phonon modes in the range 
from 50-100 meV. From the analysis of lattice dynamical 
shell model calculations, it is known that the in-plane Ni- 
0(1) oxygen stretching modes are separated in frequency 
from other vibrations, and the phonon intensity above 
65 meV is associated entirely with these vibrations. For 
the low doping samples (x — 1/8 and 1/4, not shown), 
there is little change in the ~81 meV phonon band 
but for the x — 1/3 and x = 1/2 samples, a new peak 
appears around 75 meV along with a slight hardening of 
the main band. This feature is interpreted as a splitting 
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of the 81 meV Ni-O(l) stretching band into two phonon 
bands centered at approximately 83 meV and 75 meV. 
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FIG. 1. The generalized phonon density-of-states for three 
concentrations of Laz-sSraNiCU. The frequency region shown 
(above 65 meV) consists of in-plane polarized oxygen modes 
(breathing modes) which are well-separated from other types 
of phonons. Data are offset vertically for clarity. 

In order understand this dependence of the Ni-O(l) 
stretching modes on hole concentration and its possible 
reflection of stripe ordering, we have performed a cal- 
culation of the phonon spectrum in a minimal Peierls- 
Hubbard model in 2D. Due to the strong e-1 coupling ex- 
pected in the nickelates, we resort to modeling with an in- 
homogeneous HF plus RPA numerical approach [^2] p4[ . 
This has proven to be a very robust method for study- 
ing charge localization and stripe formation, especially 
when electron-lattice coupling is strong, obviating subtle 
many-body effects and quantum fluctuations ps| . 

We use a 2D four-band extended Peierls-Hubbard 
model of a doped MO2 plane, which includes both 
electron-electron and e-1 interactions [^3|,^4). Here, for 
nickelate oxides, besides the d x 2_ y 2 orbital used in the 
cuprate oxide models [23]] , the d 3z 2_ 1 Ni d orbitals must 
be included to account for the higher spin state (S = 1) 
at half- filling (i.e. undoped). Our model Hamiltonian is 



timJn(v<ij){cl 7n( jCjn0- + H.C. 
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where c\ ma creates a hole with spin a at site i in orbital 
m (Ni d x 2_ y 2, d 3z 2_ 1 , or O p). The Ni-0 hopping t lmjn 
has two values: t p d between d x 2_ y 2 and p and ±t p d/V3 
between d 3z 2_ 1 and p. The O-site electronic energy is e p , 
and Ni-site energies are Cd and td + E z for e m , with E z 
the crystal-field splitting on the Ni site. H c describes the 
electron correlations in the Ni orbitals, 

+ J2 ( u - J / 2 ) 
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The electron-electron interactions include the on-site Ni 
Coulomb repulsions (U) as well as the Hund interaction 
(J) at the same Ni site to account for the high spin situ- 
ation. (The interplay of the two orbitals can also lead to 
pseudo Jahn- Teller distortions, but these are not our fo- 
cus here). We emphasize that, due to the large spin at the 
nickel site, Hund's rule leads to ferromagnetic exchange 
coupling —2 J, and S im = \ J] TT , cj mT CT. rT /c iTn . )T /,with er 
the Pauli matrix. For the e-1 interaction, we consider 
that the Ni-0 hopping is modified linearly by the O-ion 
displacement u i3 = u as t lmj „(uy) = f, m j„(l ± auo), 
where the + (— ) applies if the bond shrinks (stretches) 
with positive wo- For the lattice terms, we study only 
the motion of O ions along the Ni-0 bonds — other oxy- 
gen (or Ni) distortion modes can readily be included if 
necessary. It is known that for the nickelate oxides the e-1 
coupling is stronger than in cuprate oxides [pjj^ol , and 
is therefore likely to play an even more decisive role in 
the formation, localization, and nature of stripe phases. 
We adopt the following representative paremeters for the 
nickelate materials (|o): t pd = 1, A = e p — e d — 9, U = 4, 
J = 1, E z = 1 and K = 32t pd /A 2 (all in units of tpd). 
In real oxides |2fj|j2^ 1, t p d is estimated to be in the range 
1.3 eV ~ 1.5 eV. The electron-lattice coupling strength 
is varied to achieve a best fit to the neutron scattering 
data; we find a ~ 3.0. The commensurate doping cases 
are examined in a 4x4 unit supercell for x — 0, 1/2 and 
a 3x3 unit supercell for x = 1/3. Periodic boundary con- 
ditions are used. 

The densities-of-states (DOS) of in-plane phonons were 
calculated from our model at x = 0, 1/3, 1/2 and are 
shown in Fig. ^|. For the undoped the ground 

state is spatially homogeneous, only one oxygen phonon 
band appears, centered around 80.5 meV. When holes 
are added into the MO2 plane at x = 1/3, the ground 
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state is found to be a stripe pattern with more holes 
accumulating along the (1,1,0) direction, forming an an- 
tiphase domain wall within the original antiferromagnetic 
background. This is consistent with many neutron, opti- 
cal and Raman scattering experiments P,p|-[icfl . Interest- 
ingly, we find that a new phonon band appears centered 
at 75 meV. In addition, there is a hardening of the main 
phonon band (to 83 meV) corresponding to an overall 
splitting of 8 meV. From examination of the eigenvectors, 
the main character of this band is local oxygen vibrations 
in the vicinity of the stripe, i.e. having the nature of lo- 
calized "edge modes" (sec pf|] for more details). For the 
higher doping x = 1/2, the charge-ordering takes on a 
commensurate checker-board pattern. A similar split- 
ting of the 85 meV phonon band into two phonon bands 
around 83 and 75 meV is again found. At this half dop- 
ing, the checker-board ground state is in essence a com- 
mensurate charge-density-wave (CDW) system with the 
nature of an ordered binary alloy. 
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FIG. 2. The calculated densities-of-states of the oxygen 
breathing modes for various doping concentrations. The re- 
sults have been broadened with a lorentzian of width 2 meV 
for x = and 1 meV for x = 1/3 and 1/2. The electron-lattice 
coupling constant used is a = 3.0. 

The above results are in agreement with the GDOS 
data obtained from inelastic neutron scattering, which 
are shown in Fig. [l]. In addition to the splitting energy, 
even the slight hardening of the main band observed ex- 
perimentally is accounted for in the model. Besides the 
new phonon modes centered at 75 meV appearing for 
the 1/3 and 1/2 doping, another low intensity phonon 
mode around 65 meV is also predicted in our model at 
x = 1/2 (Fig. ||). The signature for these modes are 
weak, however, so that they may be difficult to detect 
in the current experiment. In so far as we have included 
only a small subset of the possible oxygen displacement 



patterns and wavevectors in the model (whereas the neu- 
tron scattering experiment samples all wavevectors and 
polarizations), the relative intensities and widths of the 
bands obtained by experiment and theory cannot be use- 
fully compared. 
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FIG. 3. The energy dependence on the electron-lattice cou- 
pling a for the 1/3 doped nickelates (solid line). For a < 2.2, 
the ground state is in an O-centered stripe phase. While for 
a > 3.0, the Ni-centered state is found as the ground state. 
The sensitive transition region is from a — 2.2 to a = 3.0. 

We emphasize that the excellent agreement between 
our model and the GDOS experimental data is achieved 
by varying the e-1 coupling strength to match the posi- 
tions of the phonon bands. As noted, the choice of a sw 3 
best fits the data; our model calculations predict a vari- 
ation of the phonon splitting with a which is quite large 
(Acu sp i it / Aa « 7 meV). Most strikingly, as illustrated 
in Fig. H an O-centered stripe is found as the ground 
state at small a < 2.2, while for a larger a > 3.0, a Ni- 
centered stripe is found as the ground state | p5| . The 
transition region includes a ~ 3.0, where the best agree- 
ment between Fig. |^ based on our model and Fig. |l| on 
the inelastic neutron scattering spectra is achieved. It 
has been suggested from various experimental data that 
stripe formation for x — 1/3 can not be simply assigned 
as Ni-centered or O-centered jl{|, but is also dependent 
on temperature. Our comparison of theory and experi- 
ment provides a possible explanation on the sensitivity of 
stripe formation; they suggest that Lai.67Sr .33NiO4 may 
be in a mixed stripe phase state, and also in a region of 
sensitivity to temperature, pressure, magnetic field, etc. 

In conclusion, we have made a study of oxygen breath- 
ing lattice vibrations in La2-^Sr ;r Ni04 via inelastic neu- 
tron scattering compared with predictions of a 2D four- 
band model, including both electron-lattice and electron- 
electron interactions. The in-plane oxygen vibrations 
above 65 meV were thoroughly investigated. The split- 
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ting of the in-plane 81 meV band upon doping into two 
subbands centered around 75 meV and 83 meV is ob- 
served experimentally and predicted theoretically, and in- 
terpreted in terms of new localized phonon modes ( "edge 
modes" at charge localized stripes). The excellent agree- 
ment between the experiment and the model strongly 
supports the view that strong electron-lattice coupling in 
this kind of material plays a decisive role on the charge 
localization and mesoscopic stripe formation. For the 
doping at x = 1/3, at which stripes are found both in ex- 
periments and our model, our results suggest that there 
may be a mixed state of O- and Ni-centered stripe phases, 
and sensitivity to temperature, pressure, and magnetic 
field. Our model calculations also predict distinctive dis- 
persion of the phonon bands, as well as inhomogeneous 
magnetoelastic coupling along the boundaries between 
charge-rich and magnetic nanophase domains |25j . These 
predictions require additional experiments for their con- 
firmation and consequences. 
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